Discontinuity of capacitance at the onset of surface superconductivity 
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The effect of the magnetic field on a capacitor with a superconducting electrode is studied within 
the Ginzburg-Landau approach. It is shown that the capacitance has a discontinuity at the onset 
of the surface superconductivity _Bc3 which is expressed as a discontinuity in the penetration depth 
of the electric field into metals. Estimates show that this discontinuity is observable with recent 
bridges for both conventional and high-Tc superconductors of the type-II. 
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Capacitors based on ferroelectric layers sandwiched be- 
tween metallic electrodes are approaching the technical 
limits of their performance. Their capacitance is not any- 
more exclusively given by the dielectric response of the 
isolating ferroelectric layer but it is reduced due to the 
penetration of the electric field into the metallic elec- 
trodes. Concerning the large scale integration of micro- 
scopic capacitors the penetration of the electrostatic field 
into the electrodes is considered as a lumped series ca- 
pacitance. From the viewpoint of fundamental research, 
however, this phenomenon offers an opportunity to study 
the interaction of metallic surfaces with an applied elec- 
tric field. 

Here we discuss a possibility to observe the penetration 
of the electrostatic field into the metal in the vicinity of 
the transition from the normal to the superconducting 
state. We focus on the third critical magnetic field Bed, 
at which field the superconducting state nucleates at the 
surface. We predict that at this field the capacitance or 
the penetration of the electrostatic field possesses a jump. 

The penetration of the electrostatic field into the nor- 
mal metal is well understood. Ku and UUman [l| have 
derived an analytic solution of the penetrating field for 
the jellium model within the Thomas-Fermi approxima- 
tion. Their simple prediction is sufficient to explain the 
experimental data [3|- Much less is known about the 
penetration of the electrostatic field into superconduc- 
tors. From the very beginning until now the history of 
this problem is full of contradictory concepts yielding a 
wide scattering of predicted values. 

The question of the penetration of the electrostatic 
field into superconductors has been firstly addressed by 
the London brothers. In their early paper in 1935 they 
have expected that the penetration depths of the electro- 
static and magnetic fields are identical [3| . One year later 
H. London measured the capacitance with superconduct- 
ing electrodes controlled by the magnetic field and con- 
cluded that the penetration of the electrostatic field into 
the metal is not changed by the transition to the super- 
conducting state [J]. While the former concept predicts 



thousands of Angstroms for conventional superconduc- 
tors, the latter concept suggests less than one Angstrom. 

Oppositly, the way from small to large penetration 
depths one meets in the papers by Anderson and cowork- 
ers. The Anderson theorem [5| states that the thcrmo- 
dynamical properties of the superconducting condensate 
do not depend on the electrostatic field. Accordingly, the 
condensate does not affect the penetration of the elec- 
trostatic field which is thus the same as in the normal 
metal. More recently, in the brief discussion of the ef- 
fect observed by Tao, Zhang, Tang and Anderson [0], the 
authors speculate about a large penetration depth of the 
electrostatic field using ideas of the Anderson model (f^ 
of the high-Tc superconductivity. 

Apparently the problem of the electrostatic field pene- 
trating the surface of the superconductor is far from being 
settled and a clear experimental message is still missing. 
In this letter we propose an experiment on the ferroelec- 
tric capacitor with one normal and one superconducting 
electrode. The magnetic field is applied to switch off the 
superconductivity and we will explore the vicinity of the 
third critical field Bcs- 

The sensitivity of ferrroelectric devices to the screening 
in metals is striking. Indeed, the typical Thomas-Fermi 
screening length in metals is about 0.5 A, while the width 
L of the insulating layer has to be about a thousand 
of Angstroms to guarantee low leakage currents. The 
direct comparison of these scales is somewhat misleading, 
however. Taking into account the dielectric constants of 
the components involved we immediately obtain 
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The ceramic ferroelectric materials have ed ^ 10"^ and 
metals have the ionic background permittivity eg ^ 4 
giving an enhancement factor ed/es ~ 250. The ca- 
pacitance can be measured with sensitivity better than 
5C /C ^ 10^^, which makes it possible to observe very 
subtle changes of the penetration depth 5L ^ lO^^A. 
First let us take a look at the interaction between the 



electrostatic field and the superconductivity. The super- 
conducting surface under the applied electrostatic field 
has been theoretically studied at various levels by Nabu- 
tovsky and Shapiro [3, M, llfl llll- They have shown 
that the phenomenological theory of Ginzburg and Lan- 
dau (GL) yields basically the same result as the mi- 
croscopic picture based on the Bogoliubov-de Gennes 
method. Their result was recovered in a simple form 
in Ref. [12], where it was shown that the effect of the 
applied electrostatic field E merely modifies the extrap- 
olation length b in the de Gennes boundary condition for 
the GL function tp, 
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This field effect is measured on the voltage scale 
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where n is the GL parameter. The logarithmic derivative 
of the critical temperature with respect to the electron 
density is of the order of unity. The need for strong 
applied fields follows from the 'relativistic' energy of an 
electron which is rather large, 111*0^ ^ 1 MeV. 

The boundary condition ([2]) restricts the solution of 
the GL equation 
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at the surface. Being non-linear, the GL equation ((4]) 
has the ability to heal any perturbation of the GL func- 
tion from its optimal value on the GL coherence length 
^ = h/ ^y2\a\m*. The boundary condition thus affects 
the GL function only in the vicinity of the surface. As 
consequence of that, the electric field has no remarkable 
effect on the bulk superconductivity. 

Saint-James and de Gennes [l3| have noticed that 
similarly to the condensation of the vapor at surfaces, 
the boundary condition ^ implies a nucleation of the 
superconducting condensate at the surface. This be- 
comes apparent at high magnetic fields, since the bulk 
superconductivity vanishes at the upper critical field Bc2 
while a thin sheet of superconducting condensate sur- 
vives near the surface up to fields Bcs ~ 1.69461 i?c2- 
Their result applies to the infinite extrapolation length, 
l/bo + E/Us = 0. 

Let us modify the method of Saint-James and de 
Gennes for a finite b or non-zero E. For the third critical 
field the GL function has an infinitesimally small ampli- 
tude so that one can neglect the cubic term in ([4]). As 
the diamagnetic current is also negligible, the vector po- 
tential reads A = {0,Bc3x,0). We have associated the 
surface with the plane x = 0. The GL equation ^ is 
then solved by the parabolic cylinder function of Whit- 
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FIG. 1: The surface critical field B^a versus the applied elec- 
tric field E (solid line). The ratio of the third to the upper 
critical magnetic field relates to the argument of the parabolic 
cylinder function, Bc^./ Bc2 = 2i/. The tangential line at zero 
bias is given as dotted line. The inset shows the field depen- 
dence of the factor in ((9|. 



taker [14 1 



,Pix,y,z)^Ne"'yD, 



2x 
T 



^kl 



(5) 



with X scaled by the magnetic length l^ — H/ieBcs) and 
ly = ~am*/he*Bc3 = P/{'2£,f- 

So far V and k are parameters of the nucleating GL 
function -0. We are looking for the solution with the 
lowest V = finin, because it corresponds to the highest 
magnetic field i?c3 for which the nucleation is possible 
at fixed temperature. This highest magnetic field Scs is 
the third critical field and its resulting value is shown in 
figure dJ 

The dependence of the third critical field Bcs on the ap- 
plied electrostatic field E indicates that the electrostatic 
field affects the surface superconductivity. The same in- 
teraction manifests itself in the effect of the magnetic 
field on the capacitance. The capacitance of the capaci- 
tor with one superconducting electrode reads 
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where S is the area of the capacitor, Cn is the capacitance 
when both electrodes are normal, and 
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is the GL free energy describing the difference between 
the normal and the superconducting states in the super- 
conducting electrode. 



Near the transition line B ^ Bcs, the GL func- 
tion has the shape of the nucleation function ([5]) with 
I' = '^min(-E'i-B) and k = k^i^{E,B). If we keep the 
amplitude A'^ as a variational parameter, the free energy 
is its biquadratic function, F^ = S{a — aE)lN'^l2/2 + 
SlN^hl{AP). Hereaij = -{he*B/m*)v^^i,,{E,B) stands 
for the kinetic energy obtained from the GL function ([5]) , 
while a = a'{T — Tc) is the temperature dependent GL 
parameter. Integrals over powers of the parabolic cylin- 
der functions are denoted as /„ — J^. dr £>"^ . The con- 
dition of minimum, dFiq/dN — 0, is satisfied by N'^^-^^ — 
[ue — a)l2/iPIi)- From F — Fn^.^^ or directly from ^ 
one obtains the free energy F = —S{aE — af'll\l{A:f5I/^. 

Now we can evaluate the jump of the capacitance, 
which appears as the magnetic field B exceeds the critical 
value -Bc3- Since aE — *• ol for B —>■ Bc3, the discontinuity 
of the inverse capacitance equals 
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where we have used daE/dE = 

— {he*B/m*){diyaiin/dE). The tangential line plot- 
ted in Fig. [T] yields dv^in/dE = -0.82 ^/C/s. The 
discontinuity in the capacitance is transparently ex- 
pressed via the discontinuity in the penetration depth of 
the electric field 
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In the rearrangement we have used 4^^ = P/iy and 
e*Bc3 = 2h/l'^. From 2v = 1.694 follows the numeri- 
cal factor 0.82'^/{2iy) = 0.397. The field dependence of 
the factor /I//4 is plotted as an inset in figured] 

To further simplify expression ^ we employ the pa- 
rameter 13 = 67r2fc|r,7(7C(3)^Fri) derived form the BCS 
theory by Gor'kov [ISJ, and rewrite it in terms of the BCS 
coherence length ^BCS = /i'yF/(l-767rfcBT'c)- Moreover we 
substitute Ug from ^ so that we obtain finally 
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We have collected all universal physical constants into 
the Bohr radius ae = 47r?i^eo/(moe^) = 0.53 A and the 
constant of fine structure e^ / {Aireohc) = 1/137. The later 
appears in the fourth power giving the very small factors 
1/137^ = 2.8x10^^. The mass of the Cooper pair is twice 
the effective mass of electrons in the metal m* — 2msmo 
and e* = 2e. For l/bo + E/Us = the factor given by the 
profile of the GL function is /I//4 = 2.42, see figure[TJ 

Equation p^ is the main result of the paper. It ex- 
presses the jump in the capacitance ([9|) in terms of mate- 
rials parameters like the logarithmic density derivative of 
the critical temperature, the coherence length ^bcs a-nd 
the GL parameter k. This result provides a convenient 



tool to access these parameters by measuring the jump 
in the capacitance at the third critical field Bcs- Indeed, 
the discontinuity is small but nevertheless observable. 

For an estimate we assume some typical numbers. The 
most sensitive measurements of capacitance performed 
in the C ~ /iF range are capable to monitor changes 
SC/C - 10-6 with error bars at SC/C - 10"^. From 
the capacitance C = eoedS/L one sees that a 1000 A- 
thick dielectric layer with ea = 10^ has an optimal area 
of 10 mm^ which is about the usual size of such samples 
[16j . The penetration depth ([9]) yields the relative change 
of the capacitance according to ([T]) . With Es — 4 and the 
above assumed values for the capacitance one finds that 
changes \SL\ > 3 x 10~6A are conveniently detectable 
with error bars oi SL ^^ 3 x lO^^A. 

It should be noted here that these estimates remain es- 
sentially valid even if the Wagner polarization diminish- 
ing effective permittivity of thin dielectric layers is taken 
into account. The expected reduction of the numbers 
above corresponds only to a factor of '--^ 2, see Ref. [iTJ. 

Now we will show that for niobium the discontinuity 
falls in the range of the error bars. For niobium at tem- 
perature T ~ 1 K one can take k '^ 1.5, see [18|, and mg — 
1.2 giving K^eJm^ = 11.7. The logarithmic derivative is 
estimated in [19| as 51nTc/i91nn = 0.74. The electron 
density n = 2.2 x lO^^/m^ yields a|n = 3.3 x 10"^ and the 
Fermi velocity wp = ?i(37r^n)^/'^/(moms) = 7.2 x 10^ m/s. 
The critical temperature T^ = 9.5 K corresponds to the 
BCS coherence length of ^bcs = 3120 A. Finally we need 
the third critical magnetic field Bed, to estimate the mag- 
netic length I. From Br.^ = 1.69 -Bc2 and the experimen- 
tal value Sc2 = 0.35 T ]^ one finds Scs = 0.59 T, which 
yields / — 325 A. With all these values we obtain from 
equation p^ the discontinuity 5L ^ 1.2 x 10^^ A, which 
is comparable to the error bar. 

There are a number of alloys [20| with the help of which 
one easily reaches a region of observable discontinuities. 
For example, 50% of niobium with 50% of tantalum has 
the critical temperature T^ = 6.25 K while the GL param- 
eter at Tc is K = 3.9 2l|. The upper critical magnetic 
field at T < Tc is 0.7T [m giving Bc3 = 1-2 T which 
yields I = 228 A. We assume that the effective mass scales 
with the GL parameter so that n/nis remains the same as 
in pure niobium along with the remaining parameters. In 
this case, the discontinuity increases to 5L = 1.1 x lO^^A, 
which is still well observable. 

We note that among intermetallic alloys there are even 
more promising candidates. The alloy of Nb-61% Ti has 
Tc = 8.95 K and k = 38.4. The upper critical magnetic 
field i?c2 = 47T corresponds to Bcs = 80 T, which is too 
high to be applied during slow measurements of the ca- 
pacitance. One has to increase the temperature for the 
measurement so that the third critical field becomes com- 
parable to a convenient field of 10 T, which corresponds 
to ; = 79A. This estimate suggests 5L = 1.5 x 10~^A 
which is fifty times larger than the experimental sensi- 



tivity. 

Detectable amplitudes of the discontinuity result also 
for high-Tc materials. Using the values of YBa2Cu307_5 
which are Tc = 90K, k = 55, to^ = 6. 92, 6= = 4, n = 
5 X lO^Vm^ ^ and dlnTc/dlnn = -2.4 ^ as well as 
I — 79 A for 10 T of the applied magnetic field, one can 
expect a discontinuity 5L = lA x 10~^A. 

It should be noted that the field effect on the high- 
Tc materials has been extensively studied within the ef- 
fort to develop superconducting devices analogous to the 
field-effect transistors [24|, |25|. There were many mea- 
surements of the field effect detecting directly changes in 
Tc with the applied electric field. These experiments em- 
ploy the largest accessible fields because the changes in 
Tc are very small. The discontinuity of the capacitance 
can supply the missing knowledge of the field effect for 
low applied fields. Since the mechanism of the field ef- 
fect on the high-Tc materials is not yet fully clarified, the 
low- field effect is of interest. 

In summary, we have shown that the capacitance of the 
planar capacitor with one normal electrode and the other 
electrode to be superconducting possesses a discontinuity 
at the third critical field Bcs- This discontinuity is large 
enough to be observed in capacitors with ferroelectric 
dielectric layers of a width of 1000 A. We would like to 
point out that compared to other regions of the magneto- 
capacitance, the discontinuity has the advantage of be- 
ing a unique feature which is not obscured by other 
properties of the insulator. Indeed, exploring strong 
electric fields one has to face the fact that the dielec- 
tric response of the ferroelectric material is non-linear. 
Scanning through temperatures one observes namely the 
Curie law of the ferroelectric transition. Moreover, the 
dielectric function of the ferroelectric isolator depends 
on the magnetic field. The measurement of the discon- 
tinuity circumvents all these problems, because all the 
troublesome dependencies are continuous at the onset of 
the surface superconductivity. 
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